Integrative spatial analysis reveals a multi-layered organization of glioblastoma.

Glioma contains malignant cells in diverse states. Here, we combine spatial transcriptomics, spatial proteomics, and computational approaches to define glioma cellular states and uncover their organization. We find three prominent modes of organization. First, gliomas are composed of small local environments, each typically enriched with one major cellular state. Second, specific pairs of states preferentially reside in proximity across multiple scales. This pairing of states is consistent across tumors. Third, these pairwise interactions collectively define a global architecture composed of five layers. Hypoxia appears to drive the layers, as it is associated with a long-range organization that includes all cancer cell states. Accordingly, tumor regions distant from any hypoxic/necrotic foci and tumors that lack hypoxia such as low-grade IDH-mutant glioma are less organized. In summary, we provide a conceptual framework for the organization of cellular states in glioma, highlighting hypoxia as a long-range tissue organizer.

Tuft cells and fibroblasts promote thymus regeneration through ILC2-mediated type 2 immune response.

The thymus is a primary lymphoid organ that is essential for the establishment of adaptive immunity through generation of immunocompetent T cells. In response to various stress signals, the thymus undergoes acute but reversible involution. However, the mechanisms governing its recovery are incompletely understood. Here, we used a dexamethasone-induced acute thymic involution mouse model to investigate how thymic hematopoietic cells (excluding T cells) contribute to thymic regeneration. scRNA-seq analysis revealed marked transcriptional and cellular changes in various thymic populations and highlighted thymus-resident innate lymphoid cells type 2 (ILC2) as a key cell type involved in the response to damage. We identified that ILC2 are activated by the alarmins IL-25 and IL-33 produced in response to tissue damage by thymic tuft cells and fibroblasts, respectively. Moreover, using mouse models deficient in either tuft cells and/or IL-33, we found that these alarmins are required for effective thymus regeneration after dexamethasone-induced damage. We also demonstrate that upon their damage-dependent activation, thymic ILC2 produce several effector molecules linked to tissue regeneration, such as amphiregulin and IL-13, which in turn promote thymic epithelial cell differentiation. Collectively, our study elucidates a previously undescribed role for thymic tuft cells and fibroblasts in thymus regeneration through activation of the type 2 immune response.

Complete human day 14 post-implantation embryo models from naive ES cells.

The ability to study human post-implantation development remains limited owing to ethical and technical challenges associated with intrauterine development after implantation1. Embryo-like models with spatially organized morphogenesis and structure of all defining embryonic and extra-embryonic tissues of the post-implantation human conceptus (that is, the embryonic disc, the bilaminar disc, the yolk sac, the chorionic sac and the surrounding trophoblast layer) remain lacking1,2. Mouse naive embryonic stem cells have recently been shown to give rise to embryonic and extra-embryonic stem cells capable of self-assembling into post-gastrulation structured stem-cell-based embryo models with spatially organized morphogenesis (called SEMs)3. Here we extend those findings to humans using only genetically unmodified human naive embryonic stem cells (cultured in human enhanced naive stem cell medium conditions)4. Such human fully integrated and complete SEMs recapitulate the organization of nearly all known lineages and compartments of post-implantation human embryos, including the epiblast, the hypoblast, the extra-embryonic mesoderm and the trophoblast layer surrounding the latter compartments. These human complete SEMs demonstrated developmental growth dynamics that resemble key hallmarks of post-implantation stage embryogenesis up to 13-14 days after fertilization (Carnegie stage 6a). These include embryonic disc and bilaminar disc formation, epiblast lumenogenesis, polarized amniogenesis, anterior-posterior symmetry breaking, primordial germ-cell specification, polarized yolk sac with visceral and parietal endoderm formation, extra-embryonic mesoderm expansion that defines a chorionic cavity and a connecting stalk, and a trophoblast-surrounding compartment demonstrating syncytium and lacunae formation. This SEM platform will probably enable the experimental investigation of previously inaccessible windows of human early post implantation up to peri-gastrulation development.

Thymic mimetic cells function beyond self-tolerance.

Development of immunocompetent T cells in the thymus is required for effective defence against all types of pathogens, including viruses, bacteria and fungi. To this end, T cells undergo a very strict educational program in the thymus, during which both non-functional and self-reactive T cell clones are eliminated by means of positive and negative selection1.Thymic epithelial cells (TECs) have an indispensable role in these processes, and previous studies have shown the notable heterogeneity of these cells2-7. Here, using multiomic analysis, we provide further insights into the functional and developmental diversity of TECs in mice, and reveal a detailed atlas of the TEC compartment according to cell transcriptional states and chromatin landscapes. Our analysis highlights unconventional TEC subsets that are similar to functionally well-defined parenchymal populations, including endocrine cells, microfold cells and myocytes. By focusing on the endocrine and microfold TEC populations, we show that endocrine TECs require Insm1 for their development and are crucial to maintaining thymus cellularity in a ghrelin-dependent manner; by contrast, microfold TECs require Spib for their development and are essential for the generation of thymic IgA+ plasma cells. Collectively, our study reveals that medullary TECs have the potential to differentiate into various types of molecularly distinct and functionally defined cells, which not only contribute to the induction of central tolerance, but also regulate the homeostasis of other thymus-resident populations.

The transcriptional and regulatory identity of erythropoietin producing cells.

Erythropoietin (Epo) is the master regulator of erythropoiesis and oxygen homeostasis. Despite its physiological importance, the molecular and genomic contexts of the cells responsible for renal Epo production remain unclear, limiting more-effective therapies for anemia. Here, we performed single-cell RNA and transposase-accessible chromatin (ATAC) sequencing of an Epo reporter mouse to molecularly identify Epo-producing cells under hypoxic conditions. Our data indicate that a distinct population of kidney stroma, which we term Norn cells, is the major source of endocrine Epo production in mice. We use these datasets to identify the markers, signaling pathways and transcriptional circuits characteristic of Norn cells. Using single-cell RNA sequencing and RNA in situ hybridization in human kidney tissues, we further provide evidence that this cell population is conserved in humans. These preliminary findings open new avenues to functionally dissect EPO gene regulation in health and disease and may serve as groundwork to improve erythropoiesis-stimulating therapies.

Exploring differential exon usage via short- and long-read RNA sequencing strategies.

Alternative splicing produces various mRNAs, and thereby various protein products, from one gene, impacting a wide range of cellular activities. However, accurate reconstruction and quantification of full-length transcripts using short-reads is limited, due to their length. Long-reads sequencing technologies may provide a solution by sequencing full-length transcripts. We explored the use of both Illumina short-reads and two long Oxford Nanopore Technology (cDNA and Direct RNA) RNA-Seq reads for detecting global differential splicing during mouse embryonic stem cell differentiation, applying several bioinformatics strategies: gene-based, isoform-based and exon-based. We detected the strongest similarity among the sequencing platforms at the gene level compared to exon-based and isoform-based. Furthermore, the exon-based strategy discovered many differential exon usage (DEU) events, mostly in a platform-dependent manner and in non-differentially expressed genes. Thus, the platforms complemented each other in the ability to detect DEUs (i.e. long-reads exhibited an advantage in detecting DEUs at the UTRs, and short-reads detected more DEUs). Exons within 20 genes, detected in one or more platforms, were here validated by PCR, including key differentiation genes, such as Mdb3 and Aplp1. We provide an important analysis resource for discovering transcriptome changes during stem cell differentiation and insights for analysing such data.

Post-gastrulation synthetic embryos generated ex utero from mouse naive ESCs.

In vitro cultured stem cells with distinct developmental capacities can contribute to embryonic or extraembryonic tissues after microinjection into pre-implantation mammalian embryos. However, whether cultured stem cells can independently give rise to entire gastrulating embryo-like structures with embryonic and extraembryonic compartments remains unknown. Here, we adapt a recently established platform for prolonged ex utero growth of natural embryos to generate mouse post-gastrulation synthetic whole embryo models (sEmbryos), with both embryonic and extraembryonic compartments, starting solely from naive ESCs. This was achieved by co-aggregating non-transduced ESCs, with naive ESCs transiently expressing Cdx2 or Gata4 to promote their priming toward trophectoderm and primitive endoderm lineages, respectively. sEmbryos adequately accomplish gastrulation, advance through key developmental milestones, and develop organ progenitors within complex extraembryonic compartments similar to E8.5 stage mouse embryos. Our findings highlight the plastic potential of naive pluripotent cells to self-organize and functionally reconstitute and model the entire mammalian embryo beyond gastrulation.

FUS-P525L Juvenile Amyotrophic Lateral Sclerosis and Intellectual Disability: Evidence for Association and Oligogenic Inheritance.

Background and Objectives: Amyotrophic lateral sclerosis (ALS) is characterized by upper and lower motor neuron degeneration, with juvenile ALS (jALS) defined as disease with age at onset (AAO) before 25 years. We aimed to identify the genetic basis of 2 unrelated patients with jALS with very rapid deterioration and early age intellectual disability (ID) and to assess association of genetic findings with both phenotypes in a large cohort of patients with ALS and controls, and in the literature. Methods: Exome sequencing was performed in 2 unrelated probands and their parents. Trio analyses included de novo, rare homozygosity, and compound heterozygosity analyses. A TaqMan genotyping assay was used to genotype ALS cohorts. A systematic literature review was conducted and additional information from authors obtained to assess prevalence of fused in sarcoma (FUS)-ALS associated with ID. Results: A de novo mutation FUS-P525L was identified in both patients. Additional variations were identified in other genes related to intellectual disabilities. Among 8 additional unrelated juvenile patients, one carried the same FUS mutation and had a similar medical history of mild ID and fulminant ALS, whereas the others did not carry any FUS coding mutations and had no reported learning or intellectual disabilities (p = 0.0083). In addition, 486 patients with ALS with AAO ≥25 years were negative for this mutation. An extensive literature review showed that among all patients with FUS-related ALS with full phenotype reports, 10.3% exhibited additional learning/intellectual disabilities. Discussion: FUS-P525L mutation was identified in 3 among 10 patients with jALS (30%) in our clinical cohort, all with a very aggressive disease course and ID. Together with literature reports, these results support a novel association between mutations in FUS and early life ID. Additional variations identified in genes related to ID and brain development in our patients (GPT2, DNAH10, and SCUBE2) may suggest a complex oligogenic inheritance for this phenotype. We propose that this mutation should be screened in patients with ALS with very early AAO, aggressive disease course, and sporadic occurrence, especially when ALS is accompanied by ID.

DestVI identifies continuums of cell types in spatial transcriptomics data.

Most spatial transcriptomics technologies are limited by their resolution, with spot sizes larger than that of a single cell. Although joint analysis with single-cell RNA sequencing can alleviate this problem, current methods are limited to assessing discrete cell types, revealing the proportion of cell types inside each spot. To identify continuous variation of the transcriptome within cells of the same type, we developed Deconvolution of Spatial Transcriptomics profiles using Variational Inference (DestVI). Using simulations, we demonstrate that DestVI outperforms existing methods for estimating gene expression for every cell type inside every spot. Applied to a study of infected lymph nodes and of a mouse tumor model, DestVI provides high-resolution, accurate spatial characterization of the cellular organization of these tissues and identifies cell-type-specific changes in gene expression between different tissue regions or between conditions. DestVI is available as part of the open-source software package scvi-tools ( https://scvi-tools.org ).

A novel mutation in TARDBP segregates with amyotrophic lateral sclerosis in a large family with early onset and fast progression.

Objective: To identify the genetic background of ALS segregating in a large Bedouin family in Israel. Methods: Exome sequencing was carried out on three siblings in a family segregating ALS, two affected and one without neurological symptoms. Filtering for causative variants and for modifiers was carried out. Eight variants were confirmed by Sanger sequencing and genotyped on nine available members of the family (three affected and six unaffected). Results: We report the identification of a novel mutation in TARDBP, p.Ala321Asp, segregating in the family. The patients are affected with early onset (average age 34.5, 21-43 years old) and fast progressive disease. The mutation is in exon 6, in the glycin-rich domain, and is predicted to be deleterious. Additional rare, potentially deleterious variants were observed in the three patients, only one of them, PLEKHG5-Phe538Leu, which is located 4.5 Mb upstream to the TARDBP, was also fully segregating in the family. Conclusion: We identified a novel mutation in TARDBP which segregates with the disease in a large family. Additional rare variants were identified, and the combination of next-generation-sequencing together with linkage analysis was optimal to identify causality and modification, emphasizing the importance of combining the two analyses. Burden of deleterious variants may be associated with early age at onset.

Rare homozygosity in amyotrophic lateral sclerosis suggests the contribution of recessive variants to disease genetics.

OBJECTIVE: to determine the occurrence of homozygous rare, in-silico damaging variants in a genetically relatively homogenous group of amyotrophic lateral sclerosis (ALS) patients. METHODS: Whole-exome-sequencing of 43 ALS patients of North-Africa Jewish origin was performed. Data were filtered to identify very rare homozygous recessive in-silico damaging variants, in genes annotated to ALS-associated cellular pathways. RESULTS: We identified a rare missense homozygous variant, p.Arg663Cys in MFN2, predicted to be damaging, in a patient with an early age at disease onset (36 years) and fast progression. An additional ALS patient carried the mutation and together established its association to ALS (p = .01). Additional homozygous variants were identified, including the risk allele p.Arg261His in NEK1, as well as variants in genes known to be associated with other neurodegenerative diseases, such as HTT (Huntington's disease), ATM (Ataxia-Telangiectasia), and ZFYVE26 (SPG15), and variants in genes previously reported as upregulated (LZTS3) or downregulated (ARMC4, CFAP54, and MTHFSD) in ALS patients. Altogether, 13 patients (30%) carried at least one homozygous rare in-silico damaging variant, of them 10 carried either another rare homozygous variant and/or a variant in a known ALS gene, which is categorized as pathogenic, likely-pathogenic or variant of uncertain significance. CONCLUSIONS: Our results suggest the contribution of recessive alleles to ALS and the possibility of burden of mutations, emphasizing the complexity of ALS genetics.

MicroRNAs and Growth Factors: An Alliance Propelling Tumor Progression.

Tumor progression requires cancer cell proliferation, migration, invasion, and attraction of blood and lymph vessels. These processes are tightly regulated by growth factors and their intracellular signaling pathways, which culminate in transcriptional programs. Hence, oncogenic mutations often capture growth factor signaling, and drugs able to intercept the underlying biochemical routes might retard cancer spread. Along with messenger RNAs, microRNAs play regulatory roles in growth factor signaling and in tumor progression. Because growth factors regulate abundance of certain microRNAs and the latter modulate the abundance of proteins necessary for growth factor signaling, the two classes of molecules form a dense web of interactions, which are dominated by a few recurring modules. We review specific examples of the alliance formed by growth factors and microRNAs and refer primarily to the epidermal growth factor (EGF) pathway. Clinical applications of the crosstalk between microRNAs and growth factors are described, including relevance to cancer therapy and to emergence of resistance to specific drugs.

EGF induces microRNAs that target suppressors of cell migration: miR-15b targets MTSS1 in breast cancer.

Growth factors promote tumor growth and metastasis. We found that epidermal growth factor (EGF) induced a set of 22 microRNAs (miRNAs) before promoting the migration of mammary cells. These miRNAs were more abundant in human breast tumors relative to the surrounding tissue, and their abundance varied among breast cancer subtypes. One of these miRNAs, miR-15b, targeted the 3' untranslated region of MTSS1 (metastasis suppressor protein 1). Although xenografts in which MTSS1 was knocked down grew more slowly in mice initially, longer-term growth was unaffected. Knocking down MTSS1 increased migration and Matrigel invasion of nontransformed mammary epithelial cells. Overexpressing MTSS1 in an invasive cell line decreased cell migration and invasiveness, decreased the formation of invadopodia and actin stress fibers, and increased the formation of cellular junctions. In tissues from breast cancer patients with the aggressive basal subtype, an inverse correlation occurred with the high expression of miRNA-15b and the low expression of MTSS1. Furthermore, low abundance of MTSS1 correlated with poor patient prognosis. Thus, growth factor-inducible miRNAs mediate mechanisms underlying the progression of cancer.

Navigator-3, a modulator of cell migration, may act as a suppressor of breast cancer progression.

Dissemination of primary tumor cells depends on migratory and invasive attributes. Here, we identify Navigator-3 (NAV3), a gene frequently mutated or deleted in human tumors, as a regulator of epithelial migration and invasion. Following induction by growth factors, NAV3 localizes to the plus ends of microtubules and enhances their polarized growth. Accordingly, NAV3 depletion trimmed microtubule growth, prolonged growth factor signaling, prevented apoptosis and enhanced random cell migration. Mathematical modeling suggested that NAV3-depleted cells acquire an advantage in terms of the way they explore their environment. In animal models, silencing NAV3 increased metastasis, whereas ectopic expression of the wild-type form, unlike expression of two, relatively unstable oncogenic mutants from human tumors, inhibited metastasis. Congruently, analyses of > 2,500 breast and lung cancer patients associated low NAV3 with shorter survival. We propose that NAV3 inhibits breast cancer progression by regulating microtubule dynamics, biasing directionally persistent rather than random migration, and inhibiting locomotion of initiated cells.

Synaptojanin 2 is a druggable mediator of metastasis and the gene is overexpressed and amplified in breast cancer.

Amplified HER2, which encodes a member of the epidermal growth factor receptor (EGFR) family, is a target of effective therapies against breast cancer. In search for similarly targetable genomic aberrations, we identified copy number gains in SYNJ2, which encodes the 5'-inositol lipid phosphatase synaptojanin 2, as well as overexpression in a small fraction of human breast tumors. Copy gain and overexpression correlated with shorter patient survival and a low abundance of the tumor suppressor microRNA miR-31. SYNJ2 promoted cell migration and invasion in culture and lung metastasis of breast tumor xenografts in mice. Knocking down SYNJ2 impaired the endocytic recycling of EGFR and the formation of cellular lamellipodia and invadopodia. Screening compound libraries identified SYNJ2-specific inhibitors that prevented cell migration but did not affect the related neural protein SYNJ1, suggesting that SYNJ2 is a potentially druggable target to block cancer cell migration.

Diurnal suppression of EGFR signalling by glucocorticoids and implications for tumour progression and treatment.

Signal transduction by receptor tyrosine kinases (RTKs) and nuclear receptors for steroid hormones is essential for body homeostasis, but the cross-talk between these receptor families is poorly understood. We observed that glucocorticoids inhibit signalling downstream of EGFR, an RTK. The underlying mechanism entails suppression of EGFR's positive feedback loops and simultaneous triggering of negative feedback loops that normally restrain EGFR. Our studies in mice reveal that the regulation of EGFR's feedback loops by glucocorticoids translates to circadian control of EGFR signalling: EGFR signals are suppressed by high glucocorticoids during the active phase (night-time in rodents), while EGFR signals are enhanced during the resting phase. Consistent with this pattern, treatment of animals bearing EGFR-driven tumours with a specific kinase inhibitor was more effective if administered during the resting phase of the day, when glucocorticoids are low. These findings support a circadian clock-based paradigm in cancer therapy.

Genome-wide mapping of IBD segments in an Ashkenazi PD cohort identifies associated haplotypes.

The recent series of large genome-wide association studies in European and Japanese cohorts established that Parkinson disease (PD) has a substantial genetic component. To further investigate the genetic landscape of PD, we performed a genome-wide scan in the largest to date Ashkenazi Jewish cohort of 1130 Parkinson patients and 2611 pooled controls. Motivated by the reduced disease allele heterogeneity and a high degree of identical-by-descent (IBD) haplotype sharing in this founder population, we conducted a haplotype association study based on mapping of shared IBD segments. We observed significant haplotype association signals at three previously implicated Parkinson loci: LRRK2 (OR = 12.05, P = 1.23 × 10(-56)), MAPT (OR = 0.62, P = 1.78 × 10(-11)) and GBA (multiple distinct haplotypes, OR > 8.28, P = 1.13 × 10(-11) and OR = 2.50, P = 1.22 × 10(-9)). In addition, we identified a novel association signal on chr2q14.3 coming from a rare haplotype (OR = 22.58, P = 1.21 × 10(-10)) and replicated it in a secondary cohort of 306 Ashkenazi PD cases and 2583 controls. Our results highlight the power of our haplotype association method, particularly useful in studies of founder populations, and reaffirm the benefits of studying complex diseases in Ashkenazi Jewish cohorts.

Tumor suppressor microRNAs: a novel non-coding alliance against cancer.

Tumor initiation and progression are the outcomes of a stepwise accumulation of genetic alterations. Among these, gene amplification and aberrant expression of oncogenic proteins, as well as deletion or inactivation of tumor suppressor genes, represent hallmark steps. Mounting evidence collected over the last few years has identified different populations of non-coding RNAs as major players in tumor suppression in almost all cancer types. Elucidating the diverse molecular mechanisms underlying the roles of non-coding RNAs in tumor progression might provide illuminating insights, potentially able to assist improved diagnosis, better staging and effective treatments of human cancers. Here we focus on several groups of tumor suppressor microRNAs, whose downregulation exerts a profound oncologic impact and might be harnessed for the benefit of cancer patients.

Context-specific microRNA analysis: identification of functional microRNAs and their mRNA targets.

MicroRNAs (miRs) function primarily as post-transcriptional negative regulators of gene expression through binding to their mRNA targets. Reliable prediction of a miR's targets is a considerable bioinformatic challenge of great importance for inferring the miR's function. Sequence-based prediction algorithms have high false-positive rates, are not in agreement, and are not biological context specific. Here we introduce CoSMic (Context-Specific MicroRNA analysis), an algorithm that combines sequence-based prediction with miR and mRNA expression data. CoSMic differs from existing methods--it identifies miRs that play active roles in the specific biological system of interest and predicts with less false positives their functional targets. We applied CoSMic to search for miRs that regulate the migratory response of human mammary cells to epidermal growth factor (EGF) stimulation. Several such miRs, whose putative targets were significantly enriched by migration processes were identified. We tested three of these miRs experimentally, and showed that they indeed affected the migratory phenotype; we also tested three negative controls. In comparison to other algorithms CoSMic indeed filters out false positives and allows improved identification of context-specific targets. CoSMic can greatly facilitate miR research in general and, in particular, advance our understanding of individual miRs' function in a specific context.